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Highly collimated, plasma- filled magnetic flux tubes are frequently observed on galactic, stellar 
and laboratory scales. We propose that a single, universal magnetohydrodynamic pumping process 
explains why such collimated, plasma-filled magnetic flux tubes are ubiquitous. Experimental evi- 
dence from carefully diagnosed laboratory simulations of astrophysical jets confirms this assertion 
and is reported here. The magnetohydrodynamic process pumps plasma into a magnetic flux tube 
and the stagnation of the resulting flow causes this flux tube to become collimated. 



The extreme collimation of astrophysical jets HI Q Q 
and the solar corona heating mechanism 4] are two seem- 
ingly unrelated astrophysical mysteries, yet both involve 
collimation of magnetic flux tubes. Astrophysical obser- 
vations 0, 0] and simulations 0, || indicate that bipo- 
lar plasma outflows (jets) are natural [3, E| features of 
young stellar objects, black holes, active galactic nuclei 
and even aspherical planetary nebula 7]. Although it 
has long been presumed [a, li| that astrophysical jets are 
magnetohydrodynamically driven, the standard models 
do not agree on a single collimation process. A similar 
issue exists in solar physics: solar spicules [lflj . promi- 
nences [Hi and coronal loops 13] are considered to 
be plasma- filled filamentary magnetic flux tubes; coro- 
nal heating models [lj, Ljjj then invoke magnetic recon- 
nection and plasma flow within such filamentary loops. 
However, the models explain neither the origin of the 
observed flows nor the extreme collimation (filamentary 
nature) of the observed structures. 

We propose that the collimation of any, initially flared, 
current-carrying magnetic flux tube is due to the fol- 
lowing process [HI: a magnetohydrodynamic (MHD) 
force resulting from the flared current profile drives ax- 
ial plasma flows along the flux tube; the flows convect 
frozen-in magnetic flux from strong magnetic field re- 
gions to weak magnetic field regions; flow stagnation then 
piles up this embedded magnetic flux, increasing the local 
magnetic field and collimating the flux tube via the pinch 
effect. Thus, the flux tube fills with ingested plasma and 
simultaneously becomes collimated. This paper presents 
direct experimental evidence for this process. We use 
ultra-high-speed imaging and Doppler measurements of 
the fast plasma flows, combined with direct density mea- 
surements before and after the filling of the flux tube. 

Our experimental setup [Tjj simulates magnetically- 
driven astrophysical jets at the laboratory scale by im- 
posing boundary conditions analogous to astrophysical 
jet boundary conditions (Fig. [lj: a disk (cathode) rep- 
resenting a central object such as a star, is coaxial and 
co-planar with an annulus (anode) representing an accre- 
tion disk. A vacuum poloidal magnetic field produced by 
an external coil links these two electrodes, mimicking a 
poloidal magnetic field threading the accretion disk. A 
radial electric field applied across the gap between the 
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FIG. 1: (log color) Typical plasma discharge sequence 
(#6577, 2 million fps, 40 ns/frame). Eight collimated 'spi- 
der legs' merge on axis to form a central column jet which 
collimates and propagates into the vacuum vessel. 



disk and annulus drives poloidal current. 

The nominal parameters of the lab experiment include 
a plasma discharge duration ~ 20 /is, poloidal magnetic 
field B ~ 0.01 - 0.3 T, gun voltage V gun - 1 - 7 kV, 
poloidal current I ~ 50 — 200 kA, density ~ io 20-22 m -3 . 
The Alfven speed is va ~ 10 5 m/s, and the ion gyroradius 
is ~ 0.5 mm, much smaller than the typical scale lengths 
(L ~ 1 — 50 cm). The plasma (3 ranges from <Cl to 
~0.5, assuming a plasma temperature ~ 10 eV, and the 
Lundquist number S = r r /rA ~ 10 3 — 10 5 , depending 
on the characteristic length used to define the resistive 
time r r (estimated from classical Spitzer resistivity) and 
the Alfven time ta- These dimensionless numbers are 
comparable to numerical MHD simulations. 

Just before breakdown, neutral gas is puffed into the 
vessel from nozzles located on the disk and the annu- 
lus. After breakdown, plasma arches form, analogous 
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FIG. 2: (color) Images from identical discharges. (a,d) 4345 
(b,e) 4346 (c,f) 4343. Nitrogen is puffed from the cathode 
nozzles and neon from the outer nozzles. Images (d-f) taken 
simultaneously with (a-c) resp., using a second camera and 
402 nm filter to image nitrogen-rich plasma. Arrow indicates 
specific spider leg analyzed in detail in Fig[7| 
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FIG. 3: Spider-leg collimation measurements, (a-c) Images 
from Figs. 0i-f. Identical white lines indicate the measure- 
ments of the spider- leg intensity, panels (d-f). Dotted line 
shows the quadratic term of the least-squares fit (solid line), 
dashed line shows the Gaussian term. FWHM (number) is 
assumed to be representative of the spider-leg width. 



to solar prominence loops [ig|. These arches are dis- 
tributed toroidally, reminiscent of spider legs, with each 
'leg' linking a gas nozzle on the disk to a corresponding 
nozzle on the annulus (Fig. [I] 4.0 [is and Fig.[2J). These 
magnetically-dominated plasma arches are initially flared 
(i.e. their minor radius increases with distance along 
their axis). Because the coil-produced poloidal vacuum 
magnetic field is weaker at the annulus than at the disk, 
the outer footpoint diameter is initially 4-5 times that 
of the inner footpoint. The eight legs quickly become 
bright and highly collimated: in 0.5 /is, as the legs fill 
up with plasma, the outer/inner footpoint diameter ratio 
has been measured to reduce by a factor of 2 (Fig. and 
can approach unity to within < 20% (Fig. |3J. The legs 
later merge (Fig. ^ 5.5-7.5 [is) to form a single central 
column (Fig. ^ 8.0 f^s), which constitutes the axially- 
expanding 'astrophysical jet' (Fig. d 8-10 [is). The cen- 
tral column can eventually kinkfFig. ^ 10-11.5 [is) and 
lead to spheromak formation 19]. The initial flaring of 
the spider legs is just the flaring of the coil- imposed vac- 
uum magnetic field, a known analytic function |2oj |. 

Measurements shown in Fig. indicate that the neu- 
tral gas density existing at breakdown is totally inad- 
equate to account for the plasma density measured in 
the spider legs or the central column. The measurements 
were taken with a custom-built fast ion gauge (FIG), with 
a 2/is response time, absolutely calibrated with a stan- 
dard commercial ion gauge. In situ measurements of the 
gas cloud in front of the electrodes were fitted numeri- 
cally to obtain a complete 3D approximation of the gas 
output from the nozzles. The numerical fit assumed a 
linear superposition of gas cones from each nozzle, each 
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FIG. 4: (color) (a, b) Measurement of flux tube flaring, on im- 
age [2f. FWHM (red marks) constant to within 20%. (c) An- 
alytical flaring of flux tube, computed from vacuum poloidal 
magnetic field 2£| . (d) Comparison of the analytical flux tube 
axis positions (white lines) with image |2t- 



with a pressure distribution having an axial exponential 
drop and a radial Gaussian profile. In particular, these 
measurements indicate that prior to breakdown, the neu- 
tral gas density at the location of a single spider leg is 
~ 10 17 m -3 (Fig. 0, which is 3-4 orders of magnitude 
less than the plasma density measured just a few [is later 
(Fig. Et). These time-resolved FIG measurements also 
show that, consistent with the molecular thermal veloc- 
ity of ~ 1 km/s, the characteristic time-scale for neutral 
density evolution is ~ 10 ms, so the neutral gas distri- 
bution is quasi-stationary during the < 50 [is plasma 
discharge duration. For a plasma density 10 3 — 10 4 times 
greater than the locally available supply of neutrals, sim- 
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FIG. 5: (color) Planar section of neutral density distribution 
just before breakdown, as measured by the fast ion gauge. 
Dashed line indicates the central axis of the experiment. 
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FIG. 6: Spectroscopic measurements: (a) center column 'as- 
trophysical jet' n e {t) inferred from Stark broadening. Solid 
line is a polynomial fit. Errors from 1-10 eV Doppler broaden- 
ing and differing Lorentzian fit parameters [2ll |. (b) Doppler 
shift of single spider-leg (N II, 399. 5nm rest-frame, dashed 
line). Solid lines are polynomial fits. Times correspond to 
end of 1 fis exposure, relative to plasma breakdown. 



pie pinching of the tube would require ~30-100 fold de- 
crease in radius. Only a 2-fold decrease is observed, so all 
spider leg particles must be ingested from the gas valves 
into the leg by MHD forces on fis time-scales. 

To visualize this fast axial flow, nitrogen was injected 
from the cathode footpoints of the spider legs and neon 
from the anode footpoints (Fig. EJ). A narrowband inter- 
ference filter (central wavelength 402 nm, Gaussian pass- 
band 10 nm) isolated the nitrogen plasma for one camera 
(Fig. EJl-f), while simultaneous, unfiltered photographs 
were obtained with a second camera (Fig. Et-c). A re- 
construction of the experimental geometry overlaid on 
the images allows comparison of the dimensions of visible 
legs with best-fit vacuum magnetic flux tube (Fig. 
The photographs indicate a brightness propagating along 
the axis of the spider- leg flux tube, from the cathode end 
to the anode end, on a 0.5 fis time-scale (Figs. ElandEJ. 
This bright front propagation occurs > 1.5 fis after the 
breakdown process has ended, demonstrating that it is 
different from streamer fronts |22| associated with break- 
down mechanisms. The time scale rules out ion acoustic 
velocities (8 km/s) from consideration and the flow direc- 
tion rules out electrostatic acceleration. The propagation 
in fact resembles a magnetized vacuum arc discharge 
and involves nitrogen plasma flow, with ~ 10 2 km/s ve- 
locities (Fig. Q) and a ~ 10 11 m/s 2 mean acceleration. 



FIG. 7: (color) Plasma flow measurements: (a) Inten- 
sity along axis (of specific spider-leg shown in Fig. EF) a t 
2.5 fis (blue diamond, Fig. EtO? 2.75 fis (green triangle, Fig. 
Et), 3 fis (red squares, Fig. EF)- Solid lines are polynomial 
fits to the data between 4 < z < 18 cm. (b) Plasma velocities 
at 2.625 fis (blue diamond line), 2.75 fis (green triangle line) 
and 2.875 fis (red square line) obtained from Eq. Q 



Flow velocities can be calculated assuming that the den- 
sity profile along the axis of the flux tube is proportional 
to the square root of the light intensity, and integrating 
the continuity equation 
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where the stagnation velocity u(L) = and L ~ 18 cm is 
the spider leg's axial length (assumed to remain constant 
during the measurements). The mean acceleration is es- 
timated from the typical velocity and time-scale (u/t ~ 
10 5 /10 -6 ~ 10 11 m/s 2 ) or from the typical length-scale 
of the flow (du z /dt ~ u 2 /(2z) ~ 10 10 /(2 x 5 x 10~ 2 ) ~ 
10 11 m/s 2 ). 

According to the model [l^, the axial component of 
the MHD pumping force is maximised on the flux tube 
axis and accelerates plasma along this axis according to 
the equation of motion 
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where a is the local flux tube radius (da/dz is the flaring), 
/ is the total current in the flux tube, u z is the axial 
plasma velocity and p is the plasma density. 

Using spider-leg flux tube dimensions measured from 
the photographs (Fig. QJ, a ~ 4 mm, da/dz ~ 0.02) 
and assuming the current inside a single leg is one 
eighth the total measured 120 kA gun current, EqEl pre- 
dicts a spider- leg plasma density of ~ 2 x 10 21 m -3 . 
This implies that within 0.5 /is, MHD forces inject 
about 2 x 10 16 nitrogen particles from the gas nozzles 
into a ~ 10 cm 3 spider-leg. Summing over all eight 
flux tubes gives a total influx pumping rate of ~ 3 x 
10 17 particles//is. Since it takes ~ 10 fis for the legs to 
merge and form the 2.5 cm diameter, 50 cm long central 
column, assuming the pumping rate is constant, the cen- 
tral column should then have a^lx 10 22 m -3 volume- 
averaged plasma density. The central column jet, located 
where 10 fis previously there was essentially vacuum, has 
ingested about 10% of the total number of particles avail- 
able from the gas feed lines. At this rate, it would take a 
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FIG. 8: Measured line spectra with Lorentzian fits: (a) single 
spider- leg, H 486.133 nm (dashed line). Exposure 0-1 /is w.r.t 
breakdown. Line of sight is perpendicular to the spider-leg, 
parallel to the cathode plane, (b) Fully formed central column 
jet (N II, dashed line is 424.178 nm rest-frame wavelength). 

few hundred /is to empty all the gas from the feed lines. 
This time is much shorter than the time-scale of fusion 
experiments, where similar filamentary flux tubes have 
been observed to emerge |2J| explosively from the hot 
plasma core and connect to the cold walls. 

The prediction for very high plasma densities in the 
spider legs and the central column has been substanti- 
ated using independent spectroscopic line measurements 
(Figs. EJl, |S|). These measurements can exhibit large 
Stark broadening indicative of high electron densities. 
For hydrogen plasmas, we measure ~ 2 x 10 20 m -3 in 
the spider-leg (Fig. [Ht) and ~ 2 x 10 22 m~ 3 in the fully- 
formed central column jet. In nitrogen discharges, central 
column jet densities 8 x 10 22 m -3 are measured (Fig. HJi 
at 4 /is, Fig. (Hb) however spider-leg densities cannot be 
resolved. Spectral lines were sampled by optical fibres 
viewing selected portions of the discharge and recorded 
with an optical multichannel analyzer. The FWHM g 
[nm] is related to n e [m -3 ] by g = 1.6 x 10 _24 n e for Nil 
H3 and g = 2 x 10- 15 n 2 e /S for H 0. The measurements 
are corrected for the instrument function and a 1-10 eV 
thermal broadening, estimated from ion line ratios. 

The pumping model also predicts that the largest 
plasma flow velocity in the tube is of the order of the 
Alfven velocity associated with the magnetic field gen- 
erated by the axial current. Here, the 15 kA current, 
4 x io 20-21 m -3 density and 4 mm typical radius of the 
spider-leg gives va ~ 70 — 220 km/s, encompassing the 
measured ~ 100 km/s shown in Fig. [Tj On the axis of a 
single spider- leg (Fig. 03), we observe blue shifts corre- 
sponding to ~ 20 — 40 km/s line-of-sight velocities. Using 
CCD images, we estimate the optical probe to be 8 cm 
along and 1 cm away from the spider- leg axis, with a line- 
of-sight angle of 72°±3°to the plasma flow direction. In 
the central column jet, the flow velocity is measured at 
~ 40 — 70 km/s (Fig. [Sb), with a line-of-sight coincident 
to the jet axis. These velocities are consistent with the 
Alfven velocity expected from the measured total 150 kA 
current and 7 x 10 22 m -3 peak density in the central col- 
umn. Both the collimation of the spider legs and the 
central column jet arise from the same process. 

The results thus provide strong evidence for flow- 
driven collimation of flared magnetic flux tubes. This 



process requires net electrical current, a supply of parti- 
cles and stagnation of the driven flow. In an astrophysical 
context, accretion disks supply plasma and flows stagnate 
at the lobes; a net axial current is understood to be nec- 
essary for self- confinement [2(j and recent observations 
of helical magnetic field structure in jets are indicative of 
net axial current [27| . The electrical circuit is then closed 
by radial currents in the accretion disk, returning via the 
cocoon. Active solar regions also exhibit unneutralized 
electrical currents 28]. We emphasize that strong flows 
and collimation can be driven by even modest amounts 
of current, before reaching the kink instability (high cur- 
rent) threshold. Supported by the US DoE and NSF. 
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